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Abstract 

We study the extension of our translationally invariant treatment of few-body nu- 
clear systems to include tensor forces and correlations. It is shown that a direct 
application of our method is not as successful for realistic V6 interactions as our 
previous results for V4 potentials suggested. We investigate the cause in detail for 
the case of 4 He, and show that a combination of our method with that of Jastrow- 
correlated wave functions seems to be a lot more powerful, thereby suggesting that 
for mildly to strongly repulsive forces such a hybrid procedure may be an appropri- 
ate description. 
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1 Introduction 



The last few years have seen an impressive progress in the microscopic descrip- 
tion of light, A < 16, nuclei directly in terms of bare realistic nucleon-nucleon 
interactions. Among the most important developments are the variational and 
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quantum Monte Carlo methods [1-4], the stochastic variational approach [5], 
the fully-converged calculations within the Faddeev and Faddeev-Yakubowsky 
theories [6], and analogous calculations in the correlated hyperspherical har- 
monics expansions [7]. 

Some of these approaches are limited by their very nature to three- or four- 
body systems. On the other hand, the quantum Monte Carlo method, the 
stochastic variational method and the correlated hyperspherical harmonics 
expansions are being extended to nuclear systems with up to around 8 or 
10 particles. At present the main limitations are both computer speed and 
computer memory, and even dealing with one additional particle requires an 
impressive growth of computing power. For example, Table 1 in Ref. [8] gives 
an indication of the computational requirements of quantum Monte Carlo 
calculations. 

At the same time we have been trying to apply what is essentially a many- 
body technique (the coupled-cluster method or CCM) to light systems. When 
dealing with finite systems of mutually interacting particles, it is important to 
separate the centre-of-mass (CM) motion exactly from the intrinsic motion. 
In the framework of the CCM it has been shown that it is possible to do so 
by using very special correlation operators [9]. Inspired by the CCM we have 
constructed a variational method that is explicitly translationally invariant, 
but which includes only a linearised correlation operator. This procedure, the 
translationally invariant configuration-interaction (TICI) method [9,10], is the 
subject of this paper. The procedure has mainly been used in a two-body 
implementation, called TICK. In recent investigations [11,12] this method 
has been extended to deal with saturated nuclei within the Op shell, using 
state-dependent (SD) pair correlations between nucleons. Calculations have 
been performed for various V4 nucleon-nucleon interactions. 

Another way to appreciate the TICK method is to consider it as a special 
formulation of the configuration-interaction (CI) method. The main difference 
with conventional CI calculations lies in the way the interacting configurations 
are selected. In the TICK approach these are restricted to the most general 
single-pair excitations which leave the CM coordinate unaffected, and are not 
limited to a few major shells, as in conventional shell model approaches. These 
configurations are very similar to those considered in the so-called potential 
basis in the hyperspherical harmonics formalism [13]. 

For systems of bosons it has been possible to extend the TICK method to 
deal with the full second-order CCM problem [14]. In this approach, called 
TICC2, the ground-state wave function contains the excitation of all possible 
independent pairs. In the specific case of 4 He with Wigner forces it was found 
that there is no significant improvement when going from TICK to TICC2, 
indicating the lack of importance of two-pair excitations by contrast with 
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triplet excitations. One should beware, however, that there is no guarantee 
that this statement holds for realistic nuclear forces with their complicated 
operatorial structure. 

In our previous work [11] we have shown how well the TICI2 method can 
deal with V4 forces and V4 correlations, defined explicitly in Sec. 2. In this 
paper we shall discuss the extension to tensor forces for Op shell nuclei. The 
TICK formalism is presented in Sec. 2, together with a hybrid method that 
adds scalar Jastrow correlations to it. In Sec. 3 we present our results, with 
a special emphasis on the case of 4 He nucleus. We find that in order to get a 
good description of the 4 He nucleus we need to mix the operatorial structure 
of the TICK wave function with a scalar Jastrow-like factor. We show that 
this improves V4 calculations as well, but not as dramatically as in the V6 
case. In Sec. 4 we draw some conclusions about the new hybrid method. 



2 Formalism 



2. 1 Linear pair correlations 

In the SUB2 truncation of the CCM, the ground-state wave function of an 
A-particle system is parametrised as 



where |<E>) is an independent-particle wave function with appropriate quantum 
numbers, and the cluster operators Si and S 2 promote one and two particles, 
respectively, out of the orbits occupied in the reference state |$) to otherwise 
empty ones. If the (uncorrelated) reference state is taken to be a Slater de- 
terminant of harmonic oscillator (HO) single-particle wave functions, the CM 
motion can be factorised exactly. The remainder of the wave function is thus 
translationally invariant. If one wishes to preserve this translational invari- 
ance in the state the excitations produced by the operators S\ and S 2 
can not be chosen independently. As has been demonstrated in Ref. [9], one 
has to construct a combination of the two operators Si and S2, denoted by 
S^' 2 \ which mixes both one- and two-body correlations in such a way that 
translational invariance is preserved. 

In the TICK method one expands the exponential of S^ 1 ' 2 **, and keeps only 
the linear term. This can be written most concisely in the coordinate repre- 



*o> =exp(S 1 + S 2 )|$), 




3 



sentation. Here the TICI2 wave function takes on the very appealing form 



A 



* (ri,...,r A ) = f( r a) $(ri,...,r A ). 



(2) 



i<j=l 



This form is very suitable for interactions that are given in coordinate space, 
such as those obtained from NN phase shifts. In order to deal with realistic 
forces we have to allow for state dependence of the pair correlation function 
/. This is done by promoting / to be an operator, and expanding it in an 
appropriate operatorial basis, /(r^) = f p {rij The basis of operators 
to be used is generally chosen to be consistent with that of the internucleon 
potential employed. For interactions with V4 structure (i.e, central, but spin- 
isospin dependent), the operators used are {6^} = {1, P CT , P T , P a P T }, where 
P a and P T are the usual spin- and isospin-exchange operators, respectively. 
When taking into account V6 interactions, we need to use a basis which also 
includes both the tensor operators Sij and SijP T . 

The uncorrected reference state is a single Slater determinant for the case of 
spin-isospin saturated systems. This can be uniquely specified by the occupied 
harmonic oscillator states. Rather than employing the standard nls basis, it 
will be convenient to work in a Cartesian basis where the states are labelled 
by three harmonic oscillator quantum numbers (n x ,n y ,n z ). The closed shell 
nuclei 4 He and 16 O are described by fully filled n < and n < 1 states; it 
also allows for the use of non-spherical reference states for the 8 Be and 12 C 



nuclei, |(0, 0, 0) 4 (0, 0, l) 4 ) and | (0, 0, 0) 4 (1, 0, 0) 4 (0, 1, 0) 4 ), respectively, which 



has proven to be better than considering a spherical reference state [11]. 

In the specific case of the 4 He nucleus, the reference wave function |<3>) is 
fully space-symmetric and so the action of the spin-exchange operator is equal 
to the isospin-exchange one, with opposite sign (i.e., space symmetry makes 
PfjP[j = — 1). Thus in this case the basis for the correlations is overcomplete 
and to avoid problems when solving for the correlation functions f p we discard 
the redundant terms, and work in a basis of only three operators, 1, P a and 
S. 

In order to determine the ground-state energy we need to compute the expec- 
tation value of both the kinetic and the potential energy operators, as well 
as that of the unit matrix, to impose normalisation. The kinetic energy and 
norm matrices are calculated straightforwardly, and we concentrate therefore 
on the more complicated matrix element of the potential energy operator 



(v) = WErty'lE^) E /(HI*o>, 



(3) 



i<j k<l m<n 
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Fig. 1. A graphical representation of the diagrams contributing to the potential 
energy, ordered according to, and labelled by, the number of particles iV involved. 
A solid line denotes the potential, and a dashed line denotes a correlation function, 
one from the bra and another from the ket state. The vertices represent particles; 
each vertex ranges over all particles in the nucleus. If the correlation functions 
occupy inequivalent positions, we get two distinct contributions from each diagram. 
This leads to the total number of 16 contributions quoted in our previous work. 



where the indices in parentheses identify the particles on which the corre- 
sponding operator acts. Many of the terms in the previous expression are 
similar, and a convenient way of dealing with them is to make a diagrammatic 
analysis [10] and get the topologically distinct diagrams to expand the matrix 
elements. The topologically distinct diagrams are schematically indicated in 
Fig. 1. 

For each diagram we have to perform a summation over all occupied single- 
particle states. Since we are dealing with saturated systems, this sum is re- 
duced to the sum over occupied harmonic oscillator states times the respective 
trace of the spin and isospin operators, where we must take account of all pos- 
sible permutations of the labels in the ket state, as explained in more detail in 
[11]. Each permutation is split into a product of transpositions with space, spin 
and isospin parts, (ij) = (i s j s ) (i a j a ){i T j T ) ■ The spin and isospin permutations 
can be represented by products of spin- and isospin-exchange operators, and 
are thus incorporated into the trace over spin and isospin. The final expression 
for matrix element related to a given diagram d involving N particles is 



(*o\r(ij)V(kl)f(mn)\* )d 

A/4 

= ^E E J2 £ P ( Tr (^j®kl®mn P Pi°...Pi% P kl...Pi T N ) 
pqr if.. .1^=1 P 

x&.(l) . . .^(7V)|/^( r ..)^( rw )/-( rmn )|0 Pif (l) . ..<f> Pifi (N)). (4) 

In this expression the 0's are the spatial parts of the single-particle wave 
functions, is the statistical weight of the diagram and Pp^^p^ is the 
representation of the spin part of the permutation (Ppq_ Pi -^ is the same for 
the isospin part). 

The computation of the required traces is greatly simplified if we work in 
the direct product representation since in this representation all the exchange 
and spin matrices have the same structure, with only one non-null element 
per row and column. In this manner the storage and matrix operations can 
be performed very efficiently and, although the number of required traces is 
large, the trace calculation may be done once and stored for later use. It is 
worth noting that the number of traces required in the V6 case is much higher 
than in the V4 case due to the complicated structure of the tensor operator. 

The last step in the calculation of the matrix elements is related to the spatial 
integrals. There is a basic difference between the V4 and V6 cases; in the 
former the operatorial structure does not contribute to the spatial part of the 
calculation, but in the latter a spatial contribution to the integrals arises. This 
is due to the appearance of the spatial coordinates in the tensor operator, and 
we must therefore consider more types of integrals. In order to perform the 
integrations we expand the correlation functions in a set of Gaussians, 

f(r)=£cfexp(-&r 2 ), (5) 
k 

with a factor r 2 for tensor components. The Gaussian expansion has proven 
to give a very accurate representation of the correlation function / when both 
negative and positive values of (5 [10] are included. In practice we use a set of 
around 10 Gaussians with fixed length parameters covering a sufficiently wide 
range of distances. The integrals are computed (see [11]) by using a recurrence 
relation for all integrations apart from those involving the coordinates in the 
potential; the remaining one-dimensional integral, which involves the poten- 
tial, is then calculated numerically. For the case of 4 He we have also performed 
calculations by representing the functions f p on a grid in coordinate space, to 
check the accuracy of our results. 

Since the expectation value of the Hamiltonian, as well as the norm, are bi- 
linear in the functions f p , we can express all expectation values as matrix 
elements, as obtained by substituting the expansion (5) in expressions such 
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as (4), and isolating the contribution from a single pair of coefficients. In this 
way we can calculate the matrix elements of the Hamiltonian 7i^ n and the 
norm M m q n . In terms of these matrices one obtains for the expectation value 
of the energy, depending on the set of coefficients {c p m } (5), the formula 

p *n_/pq q 

tj + £ C M - r> *. r pq q , (P) 
Cm J\mnOn 



where summation over repeated indices must be understood both in numerator 
and denominator separately. 

Optimising with respect to cP m * we obtain the generalised eigenvalue problem 

WM = (E + EcuWZcl (7) 



where the lowest eigenvalue is an approximation for the ground-state energy. 
We must subtract the centre-of-mass energy, 

„ ft 3 o . _ . 

E C m = — 7 « 2 , 8 
m 4 



where a is the inverse length HO parameter, from this eigenvalue to obtain the 
ground-state energy. Given that the method is strictly variational, this energy 
may also be optimised with respect to the harmonic oscillator parameter. 



2.2 Linear plus Jastrow correlations 



We shall also study a hybrid method which combines our TICI2 procedure 
with its additive state-dependent correlations, and the Jastrow method with 
its multiplicative state-independent correlations, thereby giving rise effectively 
to a correlated basis functions (CBF) description. 

It is known that Jastrow correlations are very well suited to treat the ef- 
fects of short-range repulsion. On the other hand, the TICI2 description may 
well adapt itself to represent the medium- and long-range correlations. This 
motivated us to introduce a Jastrow factor into the previous ansatz, while 
maintaining the operatorial structure of the TICI2 correlation. Incorporating 
a scalar Jastrow correlation between all pairs of nucleons in the TICK wave 
function leads to the following J-TICI2 wave function 

<Mri,...,r,0= (n^)) (EE/ P M0L]$(ri,...,^). (9) 

\i<j ) \k<l V j 
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Due to limitations in computer power, we have restricted the calculations to 
be presented below to a simplified Jastrow correlation g{r) defined in terms 
of a single Gaussian, 

g(r) = 1 +ae~ br2 . (10) 

Our wave function now depends non-linearly on three parameters, namely the 
harmonic oscillator parameter a, and the width b and depth a of the Jastrow 
function. In addition, it depends linearly on the functions f p (r) and f p (r)*. We 
have to perform a variational search of the parameters to find the minimum 
of the expectation value Ej = (^j\H\^j)/{^j\^j) - E CM . 



3 Results and discussion 

We shall investigate results obtained with the prescriptions defined in the 
previous section for several V4 and V6 interactions, and for the nuclei 4 He, 
8 Be, 12 C and 16 0. 

We shall first look at the TICI2 method for tensor forces. Afterwards, we 
discuss J-TICI2 results for the specific case of 4 He, both for V4 and V6 in- 
ternucleon interactions. At the end of this section we analyse the interplay 
of both kinds of correlations. Preliminary versions of these results have been 
presented previously in conference contributions [12,15]. 

3. 1 TICI2 results for V6 interactions 

We have performed calculations with a wide variety of interactions: the Gogny 
potential [16], the de Tourreil-Sprung super-soft-core potential (SSC) [17], and 
the Argonne potentials AV14 [18] and AV18 [19]. In all the cases we have used 
only the V6 part of the interactions, ignoring the terms related to the spin-orbit 
and other different operators. The very different character of these potentials, 
specially at short distances, can be seen in Fig. 2. Here we include, in addition 
to the potentials mentioned above, a modified version of the Reid soft-core 
potential [20,21] (called Reid-V6), used in later calculations for 4 He only, and 
discussed further in Sec. 3.3. 

The introduction of tensor correlations in our scheme increases dramatically 
the computer time needed in the calculations, as stated in Sec. 2, and it has 
taken considerable effort to deal with the case of 16 O. In Table 1 we give 
the results for the ground-state binding energy obtained for the interactions 
mentioned above. We also list the value of the harmonic oscillator parameter 
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Fig. 2. A comparison between the V6 potentials used in this paper. The horizontal 
axis in each panel shows r (in fm), the vertical axis the potential strength (in MeV). 
Each panel represents one operatorial channel, for each of the operators used to 
expand the correlation functions, as indicated in the figure. The solid line is AV14, 
the long-dashed line is AV18. The dotted line is the SSC, and the medium dashed 
line is the Gogny. The short-long dashed line is the Reid-V6 potential. 



corresponding to the minimum. Using our method, we find that the V6 parts 
of the AV14 and AV18 potentials give a very low value for the 4 He binding 
energy. The nuclei 8 Be, 12 C and 16 are even less bound, which is not what 
we expect to see. The results obtained for the SSC and especially the Gogny 
potential exhibit a more acceptable behaviour, showing at least a monotonic 
increase in binding energy for increasing values of A. 

These results are very surprising in the light of the quality of our previous 
results for V4 forces [11], and merit further investigation. One extremely im- 
portant difference is the competition between central and tensor forces. In 
the V4 forces used in our previous work most of the binding originates in the 
central channel. For all of the V6 forces used here, most binding is obtained 
from the off-diagonal tensor force, i.e., the matrix element between the central 
and tensor-correlated channel. This matrix element competes with a strong 
repulsion in the central channel. From Table 1 and Fig. 2 we see that the ori- 
gin of our problem seems to be related to the size of the repulsive core of the 
potential. We get the most sensible result for the weak Gogny force, a slightly 
worse result for the more repulsive SSC potential, and very poor results for 



9 



Table 1 

Binding energies (in MeV) of Op-shell nuclei using TICK V6-like correlations, with 
the V6 part of the Gogny, SSC, AV14 and AV18 interactions. The parameter a is 
the optimal harmonic oscillator parameter (in fm" 1 ). 



Interaction 


Nucleus 


—E 


Q 




4 He 


27.36 


0.70 


Gognv / V6 


8 Be 


40.79 


0.58 




12 C 


73.37 


0.61 




16Q 


128.68 


0.64 




4 He 


24.12 


0.68 


SSC/V6 


8 Be 


25.98 


0.54 




12 C 


39.64 


0.54 




16Q 


63.55 


0.55 




4 He 


14.77 


0.59 


At r-i a /~\ j fi 

AV14/V6 


°Be 


9.26 


0.43 




12 C 


10.50 


0.41 




16Q 


14.97 


0.40 




4 He 


15.40 


0.61 


AV18/V6 


8 Be 


11.13 


0.47 




12 C 


14.96 


0.46 




16q 


23.76 


0.46 



the AV14 and AV18 interactions. 



3.2 J-TICI2 results for V4 interactions 

In order to shed some more light on these results we have investigated the 4 He 
system more deeply. For this nucleus we can make comparisons with other 
variational calculations and with Green function Monte Carlo results. 

We have performed calculations in the 4 He system by using the J-TICI2 wave 
function of Eq. (9). With this ansatz, as already outlined in Sec. 2, we keep 
the operatorial structure of TICK but add the flexibility of the Jastrow cor- 
relations. 

In order to investigate the origin of the different quality of TICI2 calculations 
for V4 and V6 interactions, we will start by considering the simpler case of 
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Table 2 

Comparison between different results for the 4 He binding energy (in MeV) using 
several V4 potentials. S3w is the Wigner part of the S3 potential. Jastrow results 
refer to calculations using state-dependent correlations containing two Gaussians. 





TICK [11] 


J-TICI2 


SD-Jastrow [26] 


DMC [27] 


Bl 


37.86 


38.28 


38.27 [27] 


38.32 ± 0.01 


S3 


28.19 


30.16 


29.94 




S3\v 


25.41 


27.20 


27.21 


27.35 ± 0.02 


MS3 


27.99 


29.97 


29.70 




MT V 


29.45 


31.21 


30.88 


31.32 ± 0.02 


MT I/III 


30.81 


32.70 


32.01 





the V4 potentials that were studied in our previous paper [11], namely the 
Brink-Boeker Bl [22], the Afnan-Tang S3 [23] and its modified version MS3 
[24], and the Malfliet-Tjon MT I/III and MT V [25] potentials. In Table 2, we 
compare our J-TICI2 results for V4 interactions with the TICI2 ones [11], and 
with those from a recent Jastrow calculation [26] using SD correlations. In the 
case of the purely Wigner-like interactions, we also quote some diffusion Monte 
Carlo (DMC) results [27] which provide a benchmark for the calculations. 

The main difference between the J-TICI2 and the SD-Jastrow calculations 
lies in the way the state dependence of the correlations is parametrised. 
In the case of SD-Jastrow, the correlation factor is assumed to have the 
same operatorial structure as the potential, and the product of correlated 
pairs n«<j (Z) P fi ,p ( r ii)@fj) must be appropriately symmetrised to maintain the 
Fermi statistics. In the SD-Jastrow results quoted in Table 2, each of the func- 
tions g p (r) was parametrised as the sum of two Gaussians. On the other hand, 
in the J-TICI2 calculation the Jastrow correlation factor is state-independent, 
so that all of the operatorial dependence is contained in the pair correlation 
f(r). For computational reasons, the Jastrow factor in the latter case was 
parametrised by a single Gaussian, as mentioned at the end of the previous 
section. 

In these model calculations we are trying to assess the effect of superimposing 
the scalar Jastrow correlation on top of the linear pair correlations. We see 
from Table 2 that inclusion of the scalar Jastrow factor lowers the ground- 
state energy (i.e., increases the binding energy) of the pure TICI2 results by 
around 1-2 MeV. The binding energies obtained for the Wigner-like potentials 
Bl and S3 W from the J-TICI2 method and from the use of SD-Jastrow corre- 
lations are practically the same, but the former works better for more realistic 
interactions with operatorial structure. Especially interesting is the case of the 
MT potentials: for the MT V form, the J-TICI2 value is already better than 
the SD-Jastrow one, despite the fact that the former includes only one Gaus- 
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Table 3 

Same as Table 2, but referring to the V6 parts of several realistic potentials. The 
result attributed to Ref. [21] was obtained by adding the expectation values of the 
kinetic energy and the V6 part of the potential cited in that paper, and hence it is 
not a fully optimised result. 





TICK [12] 


J-TICI2 


VMC [12] 


GFMC 


Gogny/V6 


27.36 


27.58 


27.71 ± 0.06 




SSC/V6 


24.12 


26.74 


29.20 ± 0.12 




AV14/V6 


14.77 


20.37 


23.24 ± 0.08 


24.79 ± 0.20 [28] 


AV18/V6 


15.40 


21.08 


24.80 ± 0.09 




Reid-V6 


5.67 


22.70 


27.82 ± 0.12 


28.30 ± 0.12 [21] 



sian whereas the latter includes two Gaussians in the correlation functions 
g. This improvement can be explained by realizing that the MT potentials 
contain a very short-range Yukawa core. The optimal width of the Gaussian 
in the Jastrow correlation is very small (b ~ 5 fm -2 ), which cannot provide 
at the same time a good description for the long-range part of the correlation 
required to modify the independent-particle reference state. We also observe 
the proximity of the J-TICI2 and DMC results, and we hence conclude that, 
for these interactions, the J-TICI2 approach gives an excellent way to treat 
the ground state of 4 He. 

In conclusion, the V4 results indicate the relevance of the hybrid J-TICI2 
description, assigning to each part of the wave function its proper role: the 
accurate description of the short-range behaviour is assigned to the scalar 
Jastrow factor, whereas the TICK part is responsible for the medium-range 
correlations. 

3. 3 J- TICI2 results for V6 interactions 

This section deals with the most relevant application of the J-TICI2 ansatz, 
namely the consideration of realistic V6 forces. In addition to the potentials 
used before we have also considered a modified version of the original channel- 
dependent Reid soft-core (RSC) interaction [20], which presents a much harder 
core than the previous ones, and thus emphasises the role of short-range corre- 
lations. Please note that we use the same channel-independent Reid-V6 version 
of the RSC potential as has been used in the GFMC calculations of Ref. [21]. 

Table 3 summarises our V6 results. In addition to the results obtained with 
the TICK and J-TICK descriptions, the table also includes calculations with a 
V6 state-dependent Jastrow correlation computed with the variational Monte 
Carlo (VMC) method, as well as the V6 contribution to the energy com- 
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puted with the Green function Monte Carlo (GFMC) method for two interac- 
tions, AV14 [28] and Reid-V6 [21], which provide the benchmark results. Our 
state-dependent Jastrow VMC calculations have been carried out by using 
a separate Jastrow correlation function obtained by solving Schrodinger-like 
equations for each operatorial channel, as in Ref. [29]. 

The pattern of V6 results is rather different from the V4 cases. The effect of 
the inclusion of the state-independent Jastrow correlations on top of the state- 
dependent TICI2 pair correlations is most impressive for strongly repulsive 
potentials, such as the Reid-V6 [20], in which the gain in binding energy 
is about 17 MeV. Clearly, with the exception of the extremely soft Gogny 
potential, the J-TICI2 method has not converged to the level of accuracy of 
the VMC (let alone GFMC) results. We note, however, that one degree of 
freedom, which has not yet been explored due to computational limitations, is 
the use of a more general parametrisation for the Jastrow correlation function 
g in terms of multiple Gaussians. We expect some definite improvement in 
energy, but it is as yet unclear whether this will provide us with the bulk 
of the missing binding. The nature of the binding (which, as we shall see 
below, is mainly through off-diagonal components of the tensor force) is such 
that we might expect that iterations of the tensor correlations could also be 
important. The relevance of such effects would only be addressed in a full 
TICC2 calculation. 



3.4 Analysis of the results 

In this section we present our thoughts about the origins of the different be- 
haviour of the TICI2 and the J-TICI2 results. They are based on the analysis 
of the overall size of the nucleus and on a comparison of the contributions of 
the various operatorial parts to the energy. 

The first way of analysing the results is to compare the values of the harmonic 
oscillator parameter as shown in Table 4. We see that, even for V4 interactions, 
it is the increase in a, that is provided by the inclusion of the Jastrow correla- 
tion factors in Eq. (9), which drives the improvements of the energy. This need 
not surprise us; with a Jastrow correlation to kill off the short-range repulsion, 
we can use the additive CI2 correlations to get the long-range physics right. 
Without Jastrow correlations the CI2 correlation needs to play both roles, 
which works better at smaller a, where the nucleons are automatically further 
apart. The value of the HO parameter aj for the J-TICI2 ansatz does not 
change very much when considering different interactions, ranging between 
0.71 and 0.78 fm _1 and thus leads to similar overall sizes for the system. By 
contrast, the values of the Jastrow parameters present a wider range of varia- 
tion, both in depth and in range, mainly due to the fact that the short-range 
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Table 4 

TICI2 results (E, a) and J-TICI2 results (Ej, aj, a, b) for 4 He, using different 
interactions with V4 and V6 structures. Binding energies are given in MeV, the 
HO parameters a and aj in fm _1 , the depth a of the J-TICI2 Gaussian correlation 
function of Eq. (10) is dimensionless, and its width b is measured in fm~ 2 . 





TICK 




J-TICI2 






—E 


a 


-Ej 


aj 


a 


b 


V4 Interactions 














Bl 


37.86 


0.73 


38.28 


0.77 


-0.41 


1.8 


S3 


28.19 


0.72 


30.16 


0.78 


-0.70 


2.1 


S3\v 


25.41 


0.72 


27.20 


0.74 


-0.69 


2.3 


MS3 


27.99 


0.71 


29.97 


0.74 


-0.70 


2.1 


MT V 


29.45 


0.74 


31.21 


0.75 


-0.87 


5.6 


MT I/III 


30.81 


0.74 


32.70 


0.75 


-0.88 


5.0 


V6 Interactions 














Gogny/V6 


27.36 


0.70 


27.58 


0.72 


-0.29 


1.9 


SSC/V6 


24.12 


0.68 


26.74 


0.75 


-0.66 


2.1 


AV14/V6 


14.77 


0.59 


20.37 


0.71 


-0.93 


2.7 


AV18/V6 


15.40 


0.61 


21.08 


0.74 


-0.92 


3.2 


Reid-V6 


5.67 


0.43 


22.70 


0.72 


-1.05 


2.4 



behaviour of the potentials is very different from one to another. 

In Fig. 3 we compare the resulting correlation functions f p for the AV14/V6 
potential. Note that all these functions grow exponentially at large radii, but 
these contributions are moderated by the harmonic oscillator wave functions. 
For that reason we have multiplied these correlation function with the factor 
exp(— 2^-o; 2 r 2 ), arising from the relative-coordinate part of the harmonic os- 
cillator Slater determinant [11]. Sizes are related to the norm of the related 
component in the ground-state eigenfunction. As expected we see most change 
for the central correlation function. This has a node for the case of the TICK 
calculation, due to the fact that it tries to screen part of the hard-core contri- 
bution. For the J-TICI2 method this part of the correlation function is very 
smooth and structureless, apart from a small bump at small distances. We see 
that the tensor correlations are again rather smooth functions, which lends 
support to our suggestion to include only state-independent Jastrow-factors 
in our J-TICI2 method. 

A slightly deeper understanding of the mechanism of binding can be obtained 
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Fig. 3. A comparison between the TICK (solid line) and J-TICI2 (dashed line) 
correlation functions f p , multiplied with exp(— a 2 r 2 /8), for 4 He obtained from a 
solution to the generalised eigenvalue problem (7). We have chosen the V6 part 
of the AV14 potential as an example, and only plot the range for r from to 3 
fm. The normalisation of each of the correlation functions is as in the normalised 
ground-state eigenfunction, but we present separate results for each of the three 
operatorial channels. 

by decomposing the Hamiltonian in channels, where in each channel we use 
the combination of functions that was obtained from the ground-state diag- 
onalisation. This effectively reduces all matrices to three-by-three ones (one 
index for each of the central , spin-dependent and tensor-correlated channels) 
which allows for easier analysis. In a more mathematical notation, 

i<j k<l m,n 

N Pq = m E %f P *(rij) E ®lif 9 (r k iW) = E ^Zcl (11) 

i<j k<l rn,n 

where $' is either the model state $ itself or the Jastrow-correlated model 
state, $'(r 1? . . . , r^) = n«<j 9(Uj)^( r u ■ ■ ■ ■> v a), in the cases of the TICK or 
J-TICI2 calculations, respectively. The labels p and q run from one to three 
over central (1), spin (P CT ) and tensor (5), respectively. 

Since the norm matrix is diagonal, we can trivially take the square root, 
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perform the Lowdin transformation, and obtain the symmetric energy matrix, 

S = N-^HN' 1 ' 2 . (12) 



Let us compare these matrices in the case of the AV14 potential. For the 
TICK wave function the energy matrix, subtracting Ecm from the diagonal 
elements, has the value 



"-TICI2 



12.69 11.94 -96.38 
11.94 212.54 -68.09 
-96.38 -68.09 326.56 



(13) 



The eigenvector for the ground state is (0.960, 0.033, 0.278) with eigenvalue 
—14.77 MeV. Once we add the Jastrow correlation we find 



£.i- 



J-TICI2 — 



16.35 15.01 -109.66 
15.01 294.05 -74.89 
•109.66 -74.89 310.06 



(14) 



The eigenvector for the ground state is now (0.946, 0.031, 0.321) with eigen- 
value -20.37 MeV . 

In both cases the amplitude of the spin component is very small, so that 
the influence of the spin channel is not of much relevance. On the other hand, 
both diagonal matrix elements related to the central channel and to the tensor 
channel are positive, so that the central-tensor coupling is the real source of 
negative energy. 



4 Conclusions 



We have shown that the extension of our state-dependent TICK calculations 
to deal with tensor forces and tensor correlations is not as straightforward as 
our previous work suggests. The reason seems to lie mainly in the competition 
between the short-range repulsion induced by the central force and the long- 
range correlations induced by the tensor force. The presence of unscreened 
potential contributions, as shown in Fig. 1, makes the TICK method not very 
well suited to deal with realistic interactions. One might well argue that a 
full CCM calculation, where only connected diagrams appear, would alleviate 
this problem. Unfortunately, even though some progress has been made [14], 
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a coordinate-space version of the translationally invariant CCM method has 
not yet been completely and consistently formulated. 

Even if we could perform such a calculation, we would still not obtain full 
screening. Indeed, the method of Jastrow correlations is much better suited 
for such purposes. The argument about the source of repulsion (mainly of a 
central character) suggests that the J-TICI2 method proposed here should be 
an excellent starting point for describing light nuclei with realistic interactions. 
Our results for 4 He still give a few MeV less binding than the much more 
elaborate VMC calculations and the (statistically) exact GFMC calculations. 

One element that has not yet been fully explored is a better parametrisation 
for the Jastrow function, and definite conclusions about the quality of our 
approach must be delayed until we are able to deal with a more general central 
Jastrow correlation. We are also actively investigating the best way to extend 
the J-TICI2 calculations to heavier systems. 

As a last comment, we note the conceptual similarity between our approach 
and the correlated basis functions (CBF) description of many-body systems, 
which may provide another avenue for future development of hybrid methods. 
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